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Abstract 

We present a method for computing the non-perturbative mass-gap in the theory of Bosonic mem- 
branes in fiat background spacetimes. The analysis is extended to the study of membranes coupled to 
background fluxes as well. The computation of mass-gaps is carried out using a matrix regularization of 
the membrane Hamiltonians. The mass gap is shown to be naturally organized as an expansion in a 'hid- 
den' parameter, which turns out to be 4: d being the related to the dimensionality of the background 
space. We then proceed to develop a large N perturbation theory for the membrane/matrix- model 
Hamiltonians around the quantum/mass corrected effective potential. The same parameter that controls 
the perturbation theory for the mass gap is also shown to control the Hamiltonian perturbation theory 
around the effective potential. The large N perturbation theory is then translated into the language 
of quantum spin chains and the one loop spectra of various Bosonic matrix models are computed by 
applying the Bethe ansatz to the one-loop effective Hamiltonians for membranes in flat space times. 
The spin chains corresponding to the large N effective Hamiltonians for the relevant matrix models are 
generically not integrable. However, we are able to find large integrable sub-sectors for all the spin chains 
of interest. Moreover, the continuum limits of the spin chains are mapped to integrable Landau-Lifschitz 
models even if the underlying spin chains are not integrable. Apart from membranes in flat spacetimes, 
the recently proposed matrix models (hep-th/0607005 1 for non-critical membranes in plane wave type 
spacetimes are also analyzed within the paradigm of quantum spin chains. The Bosonic sectors of all the 
models proposed in (hep-th/0607005 1 are diagonalized at the one- loop level and an intriguing connection 
between the existence of supersymmetric vacua and one-loop integrability is also presented. 



1 Introduction and Summary 



In the present paper we compute the one-loop large N spectrum of various models of matrix quantum me- 
chanics describing the motion of membranes with and without supersymmetry. The purely Bosonic cases 
discussed in the paper correspond to spherical membranes in flat backgrounds which we also generalize to 
include background fluxes. We also study the one loop spectrum of the Bosonic sectors of various super- 
symmetric matrix models that were recently proposed by Kim and Park [2 1) as the regularized descriptions 
of supermembranes in non-critical super gravity theories in the background of plane wave like curved space- 
times. The matrix models corresponding to the supersymmetric membranes have explicit quadratic 'mass' 
terms in their Hamiltonians and their spectra are clearly discrete. This allows us to study their large N spec- 
trum by Bethe ansatz methods which have been very successful in the study of the the BMN matrix model. 
However, this is not the situation for the non-supersymmetric cases that we study. Since the background 
spacetime for membrane motion is chosen to be flat, the membrane Hamiltonians do not have mass terms and 
are ostensibly plagued with the problem of the existence of classical flat directions. Thus a straightforward 
application of the spin chain techniques to these models of membrane dynamics is not possible. However, the 
problem of existence of flat directions for non-supersymmetric membranes is entirely an artifact of choosing 
a particularly difficult classical theory as a starting point for quantization. The flat directions are indeed 
lifted upon quantization 23]. To make use of this simplification, we develop a method for estimating the 
dynamically generated mass terms in the quantum theory for purely bosonic membranes. The upshot of this 
technique is that one can develop a controlled large N perturbation theory both for estimating the spectral 
gap of the theories as well as for quantizing the theories around the 'quantum corrected' effective potential. 
This quantization naturally allows one to apply the techniques of quantum spin chains and the Bethe ansatz 
to systematically study the spectrum of non-supersymmetric membranes. 

The close connection between matrix quantum mechanics and theories of D branes and membranes has 
been known for quite sometime^ EJ and one of the principal motivations leading to the present work is 
the search for methods that might help us understand the gauge- theory /gravity correspondence from the 
point of view of D branes and other extended objects. A key insight gained from the recent advances in 
the AdS/CFT correspondence is that integrable structures naturally manifest themselves on both the string 
and gauge theory sides of the correspondence. As far as the spectrum of closed string like excitations, both 
from the gauge theory as well as the string theoretic point of view, are concerned, there is by now mounting 
evidence that they are described by integrable systems. On the gauge theory end, one concerns oneself with 
the spectrum of anomalous dimensions of single trace local composite operators. To the extent that it has 
been possible to check so far, this spectrum appears to be the spectrum of an integrable spin chain[3|0]. On 
the string theory side, the corresponding integrable system is nothing but the world sheet sigma model, which 
too appears to be integrable, at least at the classical level 13 E3 If one focusses on the gauge theory, then it 
is not difficult to see that the fundamental system leading to large N integrability is a quantum mechanical 
matrix model. This matrix model is nothing but the dilatation operator of the gauge theory, which in the large 
N limit can be interpreted as an integrable quantum spin chain|§l ll0llllllT2*] . Thus, keeping in mind that D 
brane and membrane world volume theories are nothing but quantum mechanical matrix models, it is perhaps 
reasonable to expect that integrable systems such as quantum spin chains should also appear naturally in 
the world volume theories of branes. Indeed, the BMN matrix model^H]) which can be interpretted as the 
theory of zero branes in type IIA theory as well as the light cone supermembrane theory of eleven dimensional 
supergravity in a plane wave background 18 is known to be integrable to rather high orders in perturbation 
theory t 20]. It provides us an example of the natural emergence of integrable systems in brane dynamics. 
However, a similar systematic approach that utilizes integrable systems has remained lacking for membrane 
motions in other backgrounds. Most notably, membrane motion in trivial/flat backgrounds have not been 
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studied from the point of view of integrable models. 

To appreciate the chief obstacle that prevents a straightforward application of quantum spin-chain tech- 
niques to matrix regularized descriptions of membranes in flat spacetimes it is worth recalling that both the 
matrix models mentioned so far, i.e the BMN model and the dilatation operator of TV = 4 SYM share a 
common feature, which is that they have discrete spectra[T51 15] . This key common feature is central to being 
able to analyze them in the large N limit as quantum spin chains. The discrete nature of the spectrum of 
the models mentioned above manifests itself very transparently at the level of the Hamiltonian in the form 
of the presence of quadratic 'mass' terms. One can then go to a basis of creation/annihilation operators and 
carry out a perturbative analysis around the oscillator vacuum. The normal ordered interaction terms can be 
understood as Hamiltonians of quantum spin chains and their contribution to the spectrum can be studied by 
applying the Bethe ansatz to the relevant spin chains order by order in perturbation theory. This connection 
between quantum spin chains and large N Hamiltonian matrix models was first worked out by Rajeev and 
Lee in[8J, and it has been extremely useful in the perturbative analysis of the spectrum and integrability 
properties of the BMN matrix model as well as the dilatation operator of M = 4 SYM |18II2(J| I3]. Presence of 
mass terms is a feature that however is not shared by various matrix models of interest. In the BMN matrix 
model m] the mass terms are a reflection of the fact that the background metric for membrane is a plane 
wave HBJ- Mat rix models for membrane motions in flat spacetimes do not have such quadratic terms. Thus 
it may seem that, at least ostensibly, membrane motions in flat spacetimes cannot be studied by the utilizing 
the connection of matrix models to quantum spin chains that has been so successful in the analysis of the 
dilatation operator of 7V=4 SYM and the BMN matrix model. This is nothing but the standard problem of 
the existence of 'flat directions' in the Hamiltonian analyses of membrane motions [TBI II 7j . However, there 
are various examples of membrane Hamiltonians whose spectra are known to be discrete at the quantum 
mechanical level even though there are apparent flat directions at the classical level. The simplest example 
of such a scenario is given by the matrix model for the motion of membranes in a d + 2 dimensional flat 
background. 



However, semiclassical analyses suggest that the spectrum of the model is indeed discrete [2"51 1251 1261 127| . In 
other words, quantum corrections lift the classical flat directions, and the quantum effective potential does 
indeed acquire a mass term. Indeed, tell-tale signs of integrable behavior in various examples of membrane 
motions in flat spacetimes have also been reported in the past literature|281 1291 13CJI Although, the 

estimation of the mass gaps of these models has been carried out using various techniques in the past, it 
would be gratifying to have a method that can enable one to do perturbative analyses around the quantum 
corrected effective potentials using the techniques of quantum spin chains. It is the development of such 
techniques that we devote the first part of the paper to. The models that we study to develop our methods 
correspond to matrix regularized membrane motions in flat background spacetimes with and without fluxes. 
The main summary of this part of the paper and the corresponding results are as follows. 

In section(3) we develop a method for computing the mass gaps in bosonic matrix models. The method 
is based on a gauge invariant rearrangement of various Feynman diagrams of the quantum mechanical model 
that allows one to compute the location of the pole of the propagator order by order in perturbation theory. 
Similar approaches for the computation of mass gaps have also been applied to various Bosonic matrix 
models as well as to supersymmetric models at finite temperature in the past, (notably in[25j). One of 
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The flat directions correspond to the configuration of commuting matrices 
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the key insights that we gain from our analysis is the emergence of a new parameter in the problem which 
controls the perturbative corrections to the dynamically generated mass term for the theory as well as 
the Hamiltonian perturbation theory around quantum corrected effective potential for the matrix models. 
Roughly speaking the parameter turns out to be 4 , where d is the number of matrices in the problem, which 
is of course related to the dimensionality of the background spacetime for the membrane theory. To put it 
differently, the obvious coupling constant of the bosonic matrix models, i.e the strength of the commutator 
squared interaction term, factors out as an overall multiplicative factor in the Hamiltonian analyses of the 
theories around the effective potential. However, one can form a dimensionless number out of the coupling 
constant of the matrix model and the dynamical mass term. It is this number that goes as 4 and appears 
to control the perturbative expansion of the matrix models. 

As mentioned previously, having a mass term in the effective potential, as well as a dynamically generated 
perturbation parameter, allows one to use the correspondence between large N matrix models and quantum 
spin chains to compute the large N spectrum of the models in perturbation theory. We carry this out in 
section(4). In this part of the paper we map the one loop effective potentials for the relevant models to 
Hamiltonians of corresponding quantum spin chains with nearest neighbor interactions. The spin chains 
corresponding to membrane motion in d + 2 dimensions have so(d) as their invariance group. This is a 
reflection of the so(d) invariance of the interaction X^,- =1 Tr[JQ, Xj] 2 term of the corresponding matrix 
models. Having a spin chain emerge as the one loop effective Hamiltonian of the matrix models is one thing, 
being able to solve it and compute the one loop spectrum of the theory is quite another. If the spin chain 
turns out to be integrable, in the sense of being part of a family of mutually commuting Hamiltonians derived 
from some underlying 1Z matrix, then one can use Bethe ansatz techniques to diagonalize the spin chains 1 . 
However, as has been explained in section(4.1) most of the so(d) invariant spin chains that emerge as the one 
loop large N Hamiltonian are not integrable in the sense described above. Nonetheless, we can still make 
progress using the following two observations. 

l:For so(2d) or so(2d + 1) spin chains with nearest neighbor interactions, there is always a rather large 
integrable sub-sector even if the spin chain is not integrable. This sub-sector corresponds the Hilbert space 
of the spin chain that contains states charged under su(d) which is contained in so(2d) or so(2d+ 1). Since 
the number of independent excitations of the spin chains is equal to the rank of the Lie algebra i.e, d — 1 for 
su(d) and d for so(2d), most of the spectrum of the spin chains of interest to us will be accessible by Bcthc 
ansatz techniques, even though the spin chains are not strictly integrable. 

2: If one were interested in the low lying excitations of the spin chain, then it is well known that they are 
captured by a classical two dimensional non-linear 'sigma' model which can be thought of as the continuum 
limit of the spin chain. For so(d) invariant spin chains we can show that the corresponding sigma model 
is always integrable even if the quantum spin chain is not. This understanding is based on realizing the 
large length (J) limit of the spin chain as a classical limit and the contraction of the quantum 1Z matrices of 
various so(d) invariant spin chains to a universal classical r matrix in the large J limit. Thus one can indeed 
understand, at least in principle, all the low lying excitations of the spin chains of interest to us irrespective 
of whether or not they are integrable. 

The final parts of chapter(4) are devoted to the one-loop perturbation theory of bosonic models with 
dynamically generated mass terms. We identify the integrable sub-sectors of the one loop spin chains and 
explicitly diagonalize them using Bethe ansatz techniques. We also give a formal derivation and proof of 
the second statement made above in section(4.4). This is carried out by deriving the non-relativistic sigma 
model for the low lying spectrum of the one loop so(d) invariant spin chains. We establish the integrability 
of the sigma models by deriving the lax pairs and monodromy matrices for them. This analysis, also brings 

1 There are several books and expository articles on Bethe ansatz techniques. For a pedagogical introduction, see|32| 
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out, in a rather explicit form, that the large J limit can be regarded as a classical limit, with -j playing the 
role of h, thus explaining the emergence of classical integrable systems in the continuum limit. 

The hnal parts of the paper are devoted to the study of supersymmetric matrix models obtained recently 
by Kim and Park from the dimensional reductions of minimally supersymmetric Yang-Mills theories in 
diverse dimensions [zj. The matrix models that we study in this section have explicit mass terms in them 
and the focus of our attention will be the large N integrability of these models. Indeed, the mechanism for 
mass generation that we utilize to study the spectra of various bosonic matrix models will not work in the 
case of most supersymmetric matrix models. The mass terms for the bosonic models arise by summing over 
certain classes of self-energy diagrams which will cancel out in the supersymmetric cases against standard 
fermionic contributions. However, there does exist a fairly large class of supersymmetric matrix models 
with explicit mass terms in their Hamiltonians; the BMN matrix model is of course an example of such 
a scenario. The BMN model is of course nothing but the light cone M(atrix) theory hamiltonian in the 
eleven dimensional plane wave background. It also has at least two other interpretations. It can also be 
thought of as the matrix regularized supermembrane Hamiltonian in the eleven dimensional plane wave 
background. As a matrix model, it can just as well also be thought of as a particular dimensional reduction 
of minimally supersymmetric Yang-Mills theory in ten dimensions (or maximally supersymmetric J\f = 4 
Yang-Mills in D=4) down to one dimension. This suggests a connection between dimensional reductions 
of supersymmetric Yang-Mill theories and super-membnrane models in one dimension higher than that of 
the original super Yang-Mills theory. It is surely instructive to probe this possible connection further. 
Minimally supersymmetric Yang-Mills theories can, other than in D = 10, be defined in dimensions six, four 
and three and two, and the matrix models obtained by the dimensional reduction of these models to one 
(time) dimension has recently been accomplished by Kim and ParkPJ. One has to be a little careful while 
relating the dimensional reductions of minimal super Yang-Mills to supermembrane theories. The naive 
dimensional reduction of super Yang-Mills to one dimension is bound to give a supermembrane theory in flat 
background. This is of course the original connection made by BFSS 13 . These matrix models do not have 
mass terms and they suffer from the problem of flat directions which cannot be cured by the method we apply 
to purely bosonic models in this paper. However, it is possible to add mass terms to the supersymmetric 
matrix models in a way that does not break any supersymmetries. The resultant mass deformed models 
can then be interpreted as supermembrane theories not in flat backgrounds, but rather in plane-wave type 
spaces. All possible supersymmetric mass deformations of the matrix models obtained from the dimensional 
reductions of minimal super Yang-Mills theories in dimensions six, four, three and two have been worked 
out by Kim and Park in their paper. The existence of the mass terms in the BMN matrix model is of course 
what makes it possible to study its large N spectrum using quantum spin chains. Moreover, the spin chains 
obtained in the pertubative expansion of the matrix model are also known to be integrable to rather high 
orders in perturbation theory. It is thus natural to ask what kind of quantum spin chains arise in the large 
N perturbation theory for the mass deformed models obtained in|21|. and whether or not these spin chains 
are integrable. This is the question that we concern ourselves with in section(5). 

Following the analysis carried out by Plefka and collaborators|2Djwe shall focus on the bosonic sectors of 
the models obtained by Kim and Park as a starting point for investigating integrable behavior in the models. 
In their paper [21], they were able to show that there are two generic types of possible mass deformations 
of the matrix models obtained from the dimensional reductions of minimal Yang-Mills theories. The mass 
deformation of the first kind corresponds to matrix models that have non-trivial maximally supersymmetric 
vacua. A second type of mass deformation, which has also been carried out by the same authors corresponds 
to supersymmetric matrix models that do not have non-trivial BPS configurations. In this paper, we are 
able to show that all the matrix models that do possess non-trivial supersymmetric vacua also correspond 
to integrable spin chains at one loop. This hints at a strong connection between integrability and super- 
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symmetry. For the mass deformed models of the second type, the one loop spin chains are not generically 
integrable. But, as mentioned previously, it is always possible to find sectors of the spin chain that do retain 
integrability, and thus give us some information about large N spectrum of the matrix models. This analysis 
is carried out in the final section of the paper. 



2 Membranes, Matrix Models and Spin Chains: A Brief Review 

In this section we gather together some known results connecting membranes to matrix models, and matrix 
models to quantum spin chains. A fuller discussion of these connections can be found in several papers. The 
relation between membrane motion and Hamiltonian matrix models is discussed in depth in |2 31 IT] while the 
use of spin chain techniques in the study of matrix models has been elaborated upon in[HI IH] • We shall refer 
to these papers for a more complete account of the results summarized in this part of the paper. 
The action of the membrane in a d + 2 dimensional flat background is given by 

S = -T J d 3 a^-deth a ^, (3) 

where T is the brane tension, which can be expressed in terms of the basic dimensional parameter in the 
problem, the Planck length, l p as 

T = 4 (4) 

h a f) = daX^dpX^ is the pull back of the spacetime metric to the membrane world volume, fj,, nu take on 
values 1 ■ ■ ■ d + 2 while a, (3 are the three world volume indices. The Nambu-Goto type membrane action can 
be replaced by a Polyakov action at the expense of the introduction of a fiducial world volume metric 7, and 
the action can be written as 

s = ¥ / (l a0 d a X^X^ - 1) . (5) 

The three diffeomorphism symmetries corresponding to the three worldvolume coordinates can be used to 
eliminate three of the six metric components, which are fixed as 

4 

7oi = 0, 700 = — ydeth. (6) 

One can introduce light front coordinates 

X ± = -L(X l ±X D+2 ), (7) 

and choose the light front gauge X + = r which eliminates two of the d+2 fields. The gauge fixed Hamiltonian 
in the light front gauge takes on the form 

H = ^ / d 2 o \X l X l + \{X\X3}{X\XJ}) . (8) 



In the above formula, X 1 are the transverse degrees of freedom and there are d of them. The curly brackets 
are the familiar Poisson brackets on the membrane surface, signifying the invariance of the action under the 
algebra of area preserving diffeomorphisms. The gauge fixed Hamiltonian can be regularized by replacing the 
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algebra of area preserving diffeomorphisms on the membrane surface by that of finite NxN matrices 23 . For 
specificity, one can take the membrane to have the topology of S 2 x R in which case one can approximate 
the algebra of area preserving diffeomorphisms on the two sphere by matrices which transform in the N 
dimensional irreducible representation of su(2). The free parameter v which was introduced in the problem 
while gauge fixing can be set to N. The details of the the approximation method leading to the regularization 
of membrane motions by the quantum mechanics of matrices can be found at various places in the literature; 
see for example 23 1 . The final answer for the regularized form of the membrane Hamiltonian in D + 2 
dimensions is given by the quantum mechanics of D hermitian matrices as 

H = .g 3 Tr(ITIT) - ^Tr ([X\ X*] [X\ X*]) (9) 

where 

g 3 = 2< (io) 

One could also allow the background spacetime to have fluxes to which the membrane degrees of freedom 
can couple. To take the simplest concrete example, one could look at zero-brane quantum mechanics of type 
IIA string theory in the presence of a non vanishing vev of the four form flux. The dynamical model is 
described by the quantum mechanics of three hermitian matrices -Xi^I], the Euclidean action for which is 

S = J dtTTr (~X? ~ l[X it X 3 ] 2 + jeykXiXjX,?) . (11) 

The cubic interaction term expresses the interaction between the matrix model and the four form flux. The 
constant four form flux is taken to be 

F$ k = -2/cey*. (12) 
T = is the zero brane tension. It is convenient to scale the matrices X; — » -4=X» and define 

9s 1 VT 
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T 

so that the Euclidean space action takes on the form 



g (is) 



S = J dtTr Qx? - tL[Xi,Xj] 2 + ^eijkXiXjX^ . (14) 

In the action the dimension of the matrices Xi is vL while those of g 2 and k are -p- and £ respectively. 
Matrix models with similar cubic couplings also arise when one considers M(atrix) theories in plane wave 
backgrounds i.e the BMN matrix model and in mass-deformed matrix models obtained from the dimensional 
reduction of super Yang-Mills theories. For the purposes of this paper, we could turn on such a four-form 
flux and couple it to the motion of a membrane in a d + 5 dimensional flat background. The flux, will pick 
out three special directions, and at the level of matrix quantum mechanics, three of the d + 3 matrices, will 
couple to the background flux. The resultant matrix mechanics will be a straightforward generalization of 
the Myers' model described above. The resultant action would be 



/ 




d^r \ - X 2 £ ^[Xi,X,-] 2 + ^ E e abc X lXj X k ] . (1 :„ 

ij—l a,b,c— 1 
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Several other examples of matrix mechanics involving cubic Chern-Simons type couplings will also be studied 
in the final parts of the paper in the context of supersymmetric matrix models. 
From Matrix Models to Spin Chains: 

Let us now briefly review the connection between the perturbative dynamics of Hamiltonian matrix models 
and quantum spin chains. This connection, which was first established by Rajeev and Lee 8 , will be central 
to the computation of the spectrum of various models of interest to us. Let us consider a matrix model 
Hamiltonian of the kind 

H = Tr Q(E? + M 2 $2) + i*£. w *^**$« + 4=M*I • jk &&$ k ) (16) 



The tensors "J 3 , \Jj 4 encode information about the cubic and quartic interaction terms. Looking ahead at the 
prospect of taking the 't Hooft large N limit in various analyses that we shall perform, we have incorporated 
various factors of -k in the generic Hamiltonian above so that it does posess a well defined 't Hooft large N 
limit. The chief qualitative difference in the nature of the Hamiltonian given above and the two examples 
discussed previously I|1UI15|I is the quadratic 'mass' term. The dynamical emergence of such mass terms is of 
course what a substantial part of the paper will be devoted to. Keeping later applications in mind we recall 
how the perturbative large N analysis of matrix model Hamiltonians becomes a problem of diagonalizing 
quantum spin chains. 

Let us start by introducing the matrix creation and annihilation operators 

Ai = -=- (pXi + iOi) (17) 

and their adjoints, in terms of which the Hamiltonian takes on the form 

H = Tr ^At% + -J^vfalAi + A^A, + A^][A k + A^][Ai + 

+ ( 2 3/ 2jV i/ V / 2) + A ^ + Ab X A * + ^)) ■ (18) 

We shall be interested in computing the first order correction to energies of single trace states such as 

\hi 2 ■■■ij>= j^Tr (AKA*> ■ ■ ■ A*') |0 > . (19) 

All such single trace states are eigenstates of the free Hamiltonian with an eigenvalue £q = /iJ. The 
perturbative corrections to the energies of the free Hamiltonian can be arranged in an expansion in powers 
of X. 

£ = /i J + -^£i + ■■■ (20) 
M 

For an eigenstate of the free Hamiltonian \T >, 

£i =<J\V\J> where V = V 4 + V 3 . (21) 

V is the vertex that contributes to the first order energy shift, and it can be built out of the cubic and 
quartic terms in (|18fl as 

V^H\V*= H ^-^ H \ (22) 
H n — J 



7 



II is the projector to the subspace of the Hilbert space orthogonal to the one spanned by states of length 
J|2(J|. while H 3 ,H 3 are the quartic and cubic terms in l|18|l respectively. 

Since only the diagonal matrix elements ofV matter for first order perturbation theory, it is sufficient 
to keep those terms in and H 3 that have an equal number of creation and annihilation operators. Not 
all such terms are important if the large N limit is taken. The terms in V that have leading order matrix 
elements of 0(1) are the ones for which normal ordering is compatible with the ordering implied by the trace, 
i.e they are terms that involve a string of creation operators followed by a string of annihilation operators 
sitting inside a trace. We can denote such operators by the symbol 0; for example, 

e;?; = iTr(At*A«,M fc ). (23) 

The action of these operators on single trace states can be written down in closed form. 

J 1 
e Ml*i*"V >= X Sl m 5l m+1 \i 1 ---i m - 1 iji m+2 ---ij > +0(jj). (24) 

m— 1 

In other words the operators 8 can be thought of as a 'machine' that runs along the entire length of the state 
and checks if any neighboring indices match the lower indices of the operators, and if they do, it replaces the 
neighboring indices by the upper indices. This leads to a map between large N matrix models and quantum 
spin chains. The single trace states can be thought of as states of a quantum spin chain with periodic 
boundary conditions, the periodicity being inherited from the cyclicity of the trace. The indices i\ ■ ■ -ij are 
to be thought of as the 'spins' of the spin chain, which take on two values. A general map between quantum 
spin chains and large N matrix models has been proposed and discussed at length in |S] and it has been 
also used in some detail to study integrability of the dilatation operator of N = 4 super Yang-Mills theory 
in|5J^]; however for the purposes of the present discussion it is sufficient to restrict ourselves to spin chains 
with nearest neighbor interactions. The so(d) invariant operators obtained from V, that survive the large 
N limit are ©H, 6}'^ andG*-'*. It is clear that O^'j acts as the number operator, thus in the large N limit 

h;;; = Tr(A^Ai) = J (25) 

The other two operators can be mapped to specific spin chain operators. The action of these spin chain 
operators on neighboring spins is 

= Pi,i+i;P{i®j)=j®i 
(->'//, = Ki,i +1 ;K(i»j) = Sij^2(m®m). (26) 

m 

In the above equations the arguments of K and P are to be interpreted as spins sitting on neighboring sites. 
P permutes neighboring spins and K traces over spin values . The large N spectrum at one loop for a 
general so{d) invariant matrix model is given by a nearest neighbor spin chain of the following generic form. 

Hi-ioop = plu+i + — o (ak,i+i + PPi.i+i + -fKi.i+i) (27) 
V I 

The coefficients a, (3, 7 carry information about the specific details about the matrix model being studied. To 
avoid various divergences that might potentially arise due to normal ordering issues, it is useful to compute 
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not the absolute energies of various states , but the energies of various states relative to a reference state. 
Since ^ h.i+i = J for all states of a given length, all terms proportional to the identity operator drop out 
when one measures the energies with respect to a particular state. The one loop spin chains that accomplish 
the task of measuring the energies with respect to a reference state are of the form 

Al = 4E^ P M+i+7^M+i) (28) 
I 

Of great interest in membrane dynamics are the following type of couplings. 

*ii«=-5(*kfy-Mi*) (29) 

so that the quartic term can be written as — &] 2 . The particular cubic coupling that will be of 
importance later in the paper is of the Chcrn-Simons type 

= iveak, (30) 

where v is a numerical constant. Keeping in mind the various models that we are going to study, we shall 

list A 1 for these two classes of models 

I: so(d) invariant matrix models for which \& = 0. 

Al = ^E(J%-%)' (3D 
II: so(3) invariant models for which ^fj k — —iveijk- 

Al ^ ^2 E - !) P M+i + \^ - 9^u+i) • (32) 



3 Dynamical Mass Generation via Resummed Perturbation The- 
ory 

We shall now revert back to the original matrix models for flat space membrane dynamics (|10I15|) and bring 
them to a form such that their spectra can be analyzed by quantum spin chain techniques. The crucial 
ingredient for this to be possible is the presence of a mass term in the matrix model Hamiltonian, the 
dynamical generation of which shall be the focus of this section. As outlined previously, we would like to 
set up a well defined perturbation theory for estimating the mass-gaps as well as the spectra of the matrix 
models of interest to us. On the face of it this might appear to be a self-contradictory goal. Dynamical 
mass generation is usually regarded as a non-perturbative phenomenon. Indeed, the methods that have 
been employed to estimate mass-gaps in matrix models in the past, see for example |25|. do indeed make use 
of non-pertrubative techniques, such as Schwinger-Dyson equations, to estimate the spectral gaps. On the 
other hand, the spin-chain techniques, that we ultimately want to make use of to do a large N Hamiltonian 
analysis appear to be an efficient method of doing the usual perturbative expansion. A way out of this 
impasse is provided is provided if one can discover a new 'hidden' parameter, which in our case will turn 
out to be , in which the spectrum of excitations as well as the mass gap of the theory can be expanded. 
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Indeed, one of the key points that we would like to bring out from this analysis is the emergence of this new 
parameter and show that it controls the perturbation theories employed both to estimate the mass gap as 
well as the one used for carrying out a Hamiltonian quantization of the theory around the quantum corrected 
effective potential. 

We start with the Euclidean action for membrane motion without fluxes H10|) . which is 

S = jdt^Tr^X i X i -^J2[X i ,X j ] 2 \.i,j = l---d-2. (33) 

We now add and subtract a mass term to the Euclidean action 

S = S m -hJ dt^m 2 Tr{X t Xi) (34) 

where 

S m = S+^ J dt^TriXiXi). (35) 

The goal is not to change the theory but only to rearrange the Feynman diagrams. Thus h will have to be 
set equal to one. We shall use the numerical value of fi only at the end of the computation. Its role is that 
of a loop counting parameter to help us carry out an efficient reorganization of the Feynman diagrams. 
The action, after appropriate insertions of factors of h reads as 

S = J dt ^ Tr + m 2 XiXi - X,j\ 2 ^j ~hj dt^TriXiXi), (36) 

Clearly, this is nothing but the original action if h is set to one. We shall let m 2 admit an expansion in h. 

m 2 = m\ + hm 2 , + h 2 m 2 + ■ ■ ■ (37) 
The propagator can be read off and it is 

»M - {(x.TMixM-9)) - >3 + J% + s(p) . 0»> 

where is the self-energy correction to the propagator. Clearly, £(p) has an expansion in powers of fi 1 

z( P ) = hz 1 ( P ) + h 2 z 2 ( P ) + --- (39) 

We can now invoke a self-consistency argument, which is that if m is indeed the gap in the spectrum i.e 
the location of the pole of the propagator then looking at the form of the propagator above, we see that 
mi , 77i2 etc must chosen to cancel the perturbative correction to the propagator that arise from the self 
energy diagrams Si, £2 ■ ■ ■ order by order in perturbation theory. This sets up a series of gap equations that 
determine the mass gap of the matrix model 2 . 

So far, we have set up the problem of determining the mass gap as a perturbative computation in fi, but 
of course, % is an invented parameter and its numerical value is one. That h is a reasonable parameter to do 

2 This generic technique is known to work rather well for gauge theories as well. For example it has been used to get a 
remarkably good estimate of the mass gap of pure Yang-Mill in 2+1 dimensions in 1241 . 
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perturbative computations in will have to be validated by an explicit demonstration that the corrections to 
the mass, when thought of as an expansion in this artificial parameter, turn out to be small in comparison 
to the leading order or one-loop estimate. That is indeed the case, and we explicitly demonstrate that in 
the next section. As a matter of fact, a two loop computation which is reported below, suggests that what 
appears to be an expansion in Ti is really an expansion in jzt- This is indeed gratifying, as this generates for 
us a bona-fide parameter in which to carry out perturbative computations of the mass gap in the problem. 
We shall now proceed to carry out explicit perturbative computations for the mass gap and arrange the 
results in a expansion. In what follows, we shall only focus on the planar diagrams, which give the 
leading order results for the mass-gap in jj ■ 



3.1 Perturbative Corrections to the Mass: 
One Loop Mass-Gap: 

The one loop correction to the mass is given by the condition that m\ cancel out the one-loop self energy 
diagram: 



fvm{ 



ftEi. (40) 



Si is given by the standard 'tadpole' integral 



w,/* ^ (41) 



2irp 2 

Ti is the standard combinatorial factor that counts the number of diagrams contributing to the tadpole 
graph in the large N limit. This number may be computed easily enough. To do that, we first observe 
that out of the two terms in the potential energy, Tr(XiXjXiXj) and Tr(XiXiXjXj), only the the second 
one contributes in the large N limit. That is the so, because, it is only in this diagram that the internal 
momentum loop of the tadpole diagram results from the Wick contraction of two neighboring X fields 
sitting inside a trace, which leads to a contribution of order N. Let us suppose, we were computing the 
correction to the < X\Xi > propagator, then clearly, the number of ways in which a X\ field from the 
Tr(XiXiXjXj)verteyi can be attached to an external line corresponding to the propagator is two. There is 
also an overall factor of d — 1 from the X fields running in the loop. Thus 

Tx = 2(d-1). (42) 

Thus the gap equation at the one loop order gives the value of the mass-gap to be 

m\ = {d-l)g z . (43) 

Two Loop Correction to The Mass: 

At the two loop level, there are two different sources of contributions to terms of 0{Ti 2 ). Since h appears 
explicitly in the formula for the propagator, the tadpole diagram, which goes as 

* I (44) 

2tt p 2 + m 2 — hm 2 

will generate a term of order h 2 . Apart from this there is the usual two loop contribution due to the 
'double-scoop' and 'sunset' diagrams 3 . Thus 

W^ - <«) 



2 TO TO 

3 We are borrowing the terminology of 42 for the Feynman diagrams 
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The first term in the expression for £2 above comes from the expansion of the formula for the 'one loop' 
tadpole diagram in higher powers of h, while the second term arises from the usual double scoop and sunset 
two loop Feynman graphs. I(p) being the total contributions of the Feynman integrals, discussed in the 
appendix. It suffices for now to know the numerical value of these diagrams, which is, 

/(p) = l((d-l) 2 + (d-l)). (46) 

(|45|l explicitly, demonstrates the intertwining of terms, that would naively be considered 'one-loop' diagrams, 
such as the 'tadpole' diagram with two loop diagrams in the second order contribution to the correction of 
the mass. Thus the two loop gap equation is: 

M^ll + h ^ d - X ) _ *V ((d _ x? + {d _ i)) = *g (47) 
m 2m 2m 4 g A 

To solve this equation perturbatively, we can set 

m 3 = g 3 {d-l)(l + hS). (48) 

The solution to 0(h) is 

5 = -2(h) m 

Thus the two loop correction to the mass can be written as 

m = ((d - l))V»,(l - ^ + • • •) = mi (l - ^ + • • •)• (50) 

Hence, we see that the two loop correction to the mass is explicitly suppressed by a factor of -jzj- Thus, as 
advertised before, the computation of the mass gap can indeed be organized in a expansion, at least to 
the two loop order. 

Even when -j^j — 1, we see that the combinatorics and the numerical values of the Feynman integrals 
involved add up to suppress the two loop contribution. As argued in |25j it is possibly too much to expect that 
the perturbation series for m actually converges, however it can possibly be regarded as a good asymptotic 
expansion, of which, we shall keep only the leading order term in later spin-chain computations. 
Aside on the Two Loop Feynman Integrals: 

The contributions of the 'double scoop' and the 'sunset' diagrams can be written as 

-4(T 2 / 2 + T 3 / 3 ) (51) 

TJl 

where T2 j3 are the combinatorial factors that count the number of diagrams contributing to the 'double 
scoop' and the 'sunset ' graphs at large N. 7 2 ,3 are the relevant Feynman integrals 

= f dpdp' 1 = I 

2 J {2TTY (p2 + l)(p'2 + l)2 8 1 ' 



and 



' ' dx f dx' — „ = — (53) 



(2^) 2 (p 2 + l)(p' 2 + l)((p-p') 2 + l) ^ 2 J-,/2 J-,/2 (tan(z) - tanOr')) 2 + 1 12 



12 



T2 is the number of Feynman graphs contributing to the double scoop diagram. This can be computed as 
follows. Once again, let us assume that we are computing the correction to the < X\X\ > propagator. Once 
again it is easy to convince oneself that the only vertex contributing to the double scoop diagram, just as in 
the case of the tadpole diagram is Tr(XiX±XjXj). The number of ways in which one of the X\ fields can 
be contracted with a an external leg of the propagator is two. The number of ways the two Xj fields can be 
contracted with the corresponding fields of the second vertex is also 2. One also picks up a factor of d — 1 
from each loop. Thus 

T 2 = 4(d-1) 2 . (54) 

In the case of the sunset diagram, we shall have to consider the cases when both the vertices are of the 
Tr(XiXiXjXj) and that when both are of the Tr(XiXjXiXj) type. Counting along similar lines as the 
case above shows that the total number of planar diagrams contributed by the first instance is 2(d — 1) and 
the number when both the vertices are of the Tr(XiXjXiXj) type is 4(d — 1). The case that involves mutual 
contractions between the vertices of these two types is of lower order in . Thus 

T 3 = 6(d-1), (55) 

leading to 

T 2 J 2 +T 3 J 3 = ~(d-l) (56) 

which was used in (|47|l . 



3.2 Mass-Gaps For Models With Chern-Simons Couplings: 

The above analysis can be easily extended to estimate the mass gap in flat space membrane motion in the 
presence of fluxes i.e to the generalization of the Meyers' model described in l|15(l . The action that we now 
consider is 

S = J dtTr Qx 2 - 9 —[X tl X 3 ] 2 + l ^e abc X a X b x)j . (57) 

In the commutator squared interaction term of the action, i,j = 1 • ■ ■ d + 3 while for the Chern-Simons 
coupling a,b,c — 1, 2, 3. The total number of Bosonic matrices in the problem is d + 3, d of which transform 
under so(d), while three of the matrices which are picked out by the anisotropy induced by the flux transform 
under so(3). The symmetry group of the problem is obviously so(d) x so(3). The one-loop mass for the 
so(d) scalars is the same as before. Adapting the result (|43|l to the present problem we can read off the one 
loop masses of the so(d) scalars to be 

m 3 so{d) =Ng 2 (d + 2). (58) 

To do a one-loop estimation of the mass gap for the so(3) scalars we shall have two Feynman diagrams 
to consider. The first one corresponds to one insertion of the quartic vertex and the second one involves 
insersions of the cubic vertex. The leading large N contribution of the insertion of the quartic vertex has 
already been discussed at length in the previous section. As far as the cubic vertices are concerned, it can 
be seen, upon carrying out the contraction of the color indices that the leading order large N contribution 
comes from the Tr{XiX2X 3 )Tr{XiX 3 X2) vertex. The momentum integrals are all convergent and can be 
evaluated in closed form by elementary means. The resulting gap equation reads as 

(d + 2)\+^—=m 3 so(3) , (59) 

m so(3) 
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where, 

A = Ng 2 (60) 

is the 't Hooft coupling of the matrix model. To solve this equation we first observe that k 2 and A 2 / 3 have 
the same dimensions. Thus, we can define a dimensionless parameter [3 such that 

k 2 = (3\ 2 /\ (61) 

In the regime where the strength of the four form flux is weak, f3 < 1, we can solve the gap equation order 
by the ansatz 

™ 3 so{3) =C\f(l3)J(p) = l + f 1 p + hl3 2 + --- (62) 

The gap equation now becomes a recursion relation for the coefficients f t while C = d + 2. The first few 
coefficients are 

1 2 

h = (d+ 2)5/3' h = 3 (( d + 2)5/3) h- ( 63 ) 

The higher coefficients in the series fall off as powers of ■ Thus the mass gaps for the so(3) scalars 
arranges itself in a k expansion, and to leading order in \ all the d + 3 scalars of the problem have the same 



4 Hamiltonian Analysis of the Matrix Models 

Having estimated the mass gaps for various Bosonic matrix models, we are now in a position to carry out a 
Hamiltonian quantization of these models and study their large N spectrum using the techniques of quantum 
spin chains. Let us begin by deriving the one-loop energy operator Ai for the most general Bosonic model 
analyzed so far, which is the so(d) x so(3) symmetric model considered in the previous section. 

After incorporating the dynamically generated mass term l|62ll . the Hamiltonian of the matrix model 
model becomes 

(1 2 2 d + 3 • 3 \ 

in 2 + HLXiXi - 9 -J2 + '-f E ^oX a X b X c . (64) 

i,j=l a.b,c—l J 

In writing down the above Hamiltonian, it is of course implied that we are starting from an effective action 
in which the self energy corrections to the propagators have been accounted for at the one-loop level. Since, 
the method of re-summed perturbation theory does not affect the vertices of the matrix model, they remain 
unchanged at the level of the Hamiltonian as well. The only modification is the inclusion of the dynamically 
generated mass term, which is the same ()(.:i2t for all the d + 3 scalars at the one loop level. 

The mass terms completely lift the classical flat directions and one can now proceed to compute the 
energy operator. As outlined in the 'introduction' the perturbative computation is best done in the basis of 
matrix creation and annihilation operators 

Ai = -=- (ml, + 01;) , A\ = -i- {mXi - iILJ . (65) 
\/2m \/2m 

The Hamiltonian becomes 

H = A^Tr (>% - —L^ + A* A 3 + A^f + i —^^ ^JU + A^){A 3 + A^)(A k + A*) 

(66) 
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where fi is defined to be 



m = A 1 / 3 M = A 1 / 3 (C/(/?)) 1 / 3 . 



(67) 



Since the symmetry of the problem splits into so(d) x so(3) it makes sense to write down the energy operator 
for the so(3) and the so(d) part of the Hamiltonian separately. That can be done easily enough using the 
basic formulae (|31I32|I : and the result is 

^o(d) =-^E - ^M+i) • (68) 

A 1)So(3 ) = 2^2 E - !) P ^+i + \^ ~ M 2 ) K u+^) * ^ = |r ( 69 ) 

It is worth noting that fi is essentially the same quantity that controls the perturbative expansion of the 
mass-gap discussed earlier. From the above formulae, it is also clear that /i is the dimensionless parameter 
that controls the Hamiltonian perturbation theory as well. 



4.1 Integrability: 

Having a spin-chain formalism for the one-loop energy operator does not automatically imply that we shall 
be able to solve for its spectrum. The obvious obstacle in this direction is that generic so(d) invariant 
quantum spin chains are not integrable by the methods of the algebraic Bethe ansatz. For an so(d) invariant 
quantum spin chain, with nearest neighbor interactions, d > 2 and spins in the defining representation, to 
be integrable, the Hamiltonian has to be of the following formal ES] 

H = aJ2(9'Pl,l+i + Ku+i) , 9' = (i - ")■ (70) 

a is an arbitrary coupling constant. One is of course free to add a term proportional to the identity operator 
to the Hamiltonian and not loose integrability. It is easy to see that neither one of I|68I69|) is of the form 
required by integrability. However, one can still make progress, as non-integrable spin chains also contain 
sectors can indeed be studied by the Bethe ansatz. To proceed further, let us briefly recall the construction 
of integrable spin chains starting from an underlying TZ matrix, which for so(d) spin chains ispfl I33j 

K 12 (u) = u(hg' - u)I u + h(hg' - u)P X2 + huK 12 . (71) 

g' = (1 - f ). For later uses, we have put in a parameter h in the 1Z matrix to act as a bookkeeping device. 
The numerical value of fi is one and it is nothing but the strength of non-commutativity of the spin operators 
at the same lattice point. The 1Z matrix is an operator that acts on a tensor product of two vector spaces 
both of which are C d when the spins are in the defining representation of so(d). The indices 1,2 refer to 
the respective vector spaces, while u is the spectral parameter 4 . For a spin chain with J sites the transfer 
matrix is built out of the 1Z matrix as 

T{u) = 1Z i(u)lZ 02 (u) ■ --IZojIu). (72) 

4 For a pedagogical introduction to the use of TZ matrices in the construction of integrable spin chains, see 1321 . 
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The subscript refers to the auxiliary space, which is not a part of the physical Hilbert space of the spin 
chain, and it will be traced over. Let Ttq denote the trace over the auxiliary space. The transfer matrix can 
be obviously be expanded in powers of the spectral parameter 



T(u)=5> m T m (73) 



and taking the trace over the auxiliary space of both sides of the above equations produces a generating 
function for the local conserved charges of the spin chain. 

t(u)=Tr T(u) = Y / u m t m . (74) 

I a 

Finding a closed form of for the local charges is usually a tedious but straightforward computation; however, 
we shall be interested in only the first few charges, especially the second one which corresponds to the 
Hamiltonian. The first charge is 

J 

to = Tr l[(h 2 g r ) J P 0l (75) 
i=i 

It can be shown that this charge corresponds to the discrete lattice shift operator. 

*o(Vi <8> V 2 <8 • • • ® Vj) = (h 2 g') J (V 2 ® V 3 ® • • • V/ ® Vi). (76) 
The next term in the expansion is 

J 

h = X> {hg'Pi,i+i - + hKi,i+i) . (77) 

1=1 

Since to is a constant of motion and we are free to add constant terms to the Hamiltonian it is easy to see 
that the most general form for the nearest neighbor Hamiltonian with a so(n) symmetry is 

J 

H = aJ2 (9'Pi,i+i ~ Ph,i+i + Ku+i) (78) 
i=i 

where a and (3 are arbitrary constants. The crucial requirement for integrability is the relative coefficient 
between the permutation and trace operators which has to be (1 — |) for an so(d) spin chain. This clearly 
indicates the lack of integrability of the energy operators (|68I69|I . 
Integrable Sub-sectors: 

Even when an so(d) spin chain is not integrable it is possible to find integrable sub-sectors of the theory. As 
has been reviewed in the appendix, if one constructed and su(d') integrable spin chain out of the relevant 1Z 
matrix the integrable nearest neighbor Hamiltonian would be 



# = a J^ify+i. (79) 



i 



so(d = 2d') or so(d = 2d' + 1) of course contains su(d'). If we focussed only on spins transforming under 
this su(d'), then at the level of the Hamiltonian, it would correspond to dropping the trace operator, and 
we shall be left with an integrable su(d') sector of the theory. To find this embedding of su(d') at the level 
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of the Hilbert space, one needs to isolate a subspace of the Hilbert space of the spin chain such that 
1: It is annihilated by the trace operator 

2: It remains closed under the action of the spin chain Hamiltonian. 

When d — 2d', one can always find an integrable su(d!) sub-sector for the theory as follows. One can form 
d' complex combinations of the 2d' spins i% - • • iid 1 as 

X x = i x + -s/— Tid'+i = id' + y/—U2d>- (80) 

The sub-space of the spin-chain Hilbert space spanned by states such as 

\f>=f h ... lj \l h ---l h > (81) 

involve only traceless so(d) tensors and are hence annihilated by the trace operator. In this sub-space, the 
spin chain Hamiltonian involves only the identity and the permutation operators and is hence integrable. 
When d = 2d' + 1, the biggest closed integrable sector is the same as the case when d = 2d' . Of course, 
we cannot use this simplification to diagonalize any part of the Hamiltonian 1)69(1 . for which we shall resort 
to some approximate methods described later in the paper. However, the presence of the integrable su(d') 
sector allows us to study the spectrum of Bosonic membranes in flat backgrounds in dimensions greater than 
four. 

There is an important exception to the list of integrable spin chains formed out of so(d) invariant 1Z 
matrices discussed above and that corresponds to d = 2. In four dimensions, which corresponds to the case 
d = 2, the above analysis of integrability does not apply and the corresponding one-loop spin chain does turn 
out to be integrable! Since the Bethe ansatz for this particular case is rather transparent and the generic 
nature of the spectrum of the d — 2 case is indeed the same as that of the spin chains corresponding to 
integrable su(d') subsectors of membranes in higher dimensions we shall start by analyzing this special case 
first. We shall then generalize the result to the integrable su(d') sectors of various bosonic membranes at 
one loop. 



4.2 The Special Case of d = 2 

To solve for the spectrum of the special cased d = 2, the one loop energy operator 1(68(1 we shall employ the 
techniques of coordinate space Bethe ansatz. Let us first write the trace and permutation operators in terms 
of the Pauli spin matrices as follows. 



l - ^1(1) (8 1(1 + 1) + ® + 1)J 



n,i+i 2 

K u+ i = ^{I(l)®I{l + l) + a x {l)®a%l + l)+a z (l)®a z {l + l)-cjy{l)®ay{l + l)). (82) 
The energy operator can then we expressed as 

A Mo(2) = ± (jy*{l)<7*(l + 1) + o*{l)o*{l + 1) + 3ay(l)ay(l + 1)]^ (83) 

This is nothing but the anisotropic Heiseberg chain, also known as the XXZ spin chain 5 . The spin chain has 
an obvious symmetry which allows one to carry out rotations in the x — z plane by rotating a x and a z among 

5 The name the model refers to the fact that the direction of the anisotrpy is usually chosen to be the Z direction. We can 
bring the model to that form simply by re-labeling the sigma matrices. 
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each other. This is nothing but a reflection of the original so(2) symmetry of the interaction Tr[Xi,X% 
term of the matrix model Hamiltonian. We can also express the hamiltonian in the form 



A Mo(2) = 7 (j2[Sh(l)S%(l + 1) + S?(0S 2 x (i + 1) + S(Sl(l)S 2 2 (l + 1) + Sf(l)Sj(l + 1))]^ (84) 

In the above equation, 7 = and 5 = —3. We have introduced the Weyl operators, Sj(l), that check if 

the spin at the I th site is equal to j and if it is, then Si (I) replaces the value of the spin by i. If the value 
of the spin at the I th site does not match with j then this operator annihilates the state. These operators 
satisfy the associative Weyl agebra 

Sj(Z)S£(0= $Si(l) (85) 

at the same site and the Weyl operators at different sites commute with each other. With this form (|84|1 of 
the Hamiltonian, it is clear that the Ferromagnetic state, with all the spins up(l) is a ground state of the 
Hamiltonian. 

|0>=|1,1, !•••!> (86) 

zero energy. The complete set of excited states of the Hamiltonian can be built by considering linear 
combinations of states with a fixed number of down/impurity spins. The one magnon state or the state with 
a single impurity can be taken to be of the form 



>= £/(0|Z >, (87) 



1=1 

where \l > is the state with an upturned spin at the Ith site. If the state is taken to be a plane wave 

f(l) = e** (88) 

then it is easy to see that it will be a eigenstate of the Hamiltonian with eigenvalue 

£(p) = -4 7 (l + 2sin 2 (|)). (89) 

The allowed values of the momenta are determined from the condition that the state be cyclically symmetric 
i.e 

e^ J = l^p=^. (90) 

The one magnon state gives us the dispersion relation l|9(J|) and we shall have to consider the two magnon 
state to get the two magnon scattering matrix. The state with two 'down' spins can be written as 

J 

1*2 >= f(xi,X 2 )\x U X 2 > (91) 

(xi<x 2 )=l 

where it is understood that \xi, X2 > is the state with upturned spins at xi and X2- Once again we can make 
a Bethe ansatz by taking a linear combination of plane waves that interact with each other by exchanging 
momenta 

f(xi,x 2 ) = e l{piXl+P2X2) + S{p llP 2)e l ( p2Xl+PlX2 \ (92) 
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Such a state will be an eigenstate of the Hamiltonian if the two magnon scattering matrix S(pi,p 2 ) is given 
by 

I + e i(pi+ P2 ) + 25e^ 
S{ P^) = - i + e*(pi+P2) +2 fc'Pi - (93) 

The eigenvalue is given by 

E(pi,p 2 ) = -4 7 ]T (l + sin2 (f )) > ( 94 ) 

i=l,2 

which shows that the energies of the two magnons simply add up. The momenta are determined by the 
condition that the total phase picked up by one of the magnons by travelling around the spin chain is nothing 
but the two body scattering matrix, i.e. 

e ipiJ = S(pi,p 2 ),e ip2j = 5(p 2 ,Pi). (95) 

These are the Bethe equations for the two unknowns pi , p 2 and they provide us with the complete solution 
to the two magnon problem. The XXZ quantum spin chain is an integrable system, and one of the 
manifestations of integrability is the factorized nature of the multi-magnon scattering matrix into two particle 
S matrices. In other words, the above equations also provide us with the solution to the the m < J magnon 
problem. The eigenstate with m down spins has energy given by 



£ = -4 7 £(l + sin 2 (f)), (96) 
i=i 

with the momenta determined by 

m 

e^ J = J] S(p k , Pj ). (97) 

These are m equations for the m unknown momenta that provide us with the complete spectrum of the spin 
chain. To get a feel for the qualitative nature of the spectrum as a function of J we can invoke a BMN 
limit and work in a the dilute gas approximation, where the number of magnons is very small compared to 
the length of the chain i.c m « J. We can write l|97() by taking the logarithm of both sides of the equation 
as 

p k J = 2mr + G(jp 1 ,p 2 ). (98) 

is the two body scattering phase shift. In the dilute gas approximation, the magnons do not scatter at 
all, and have 0^0. The dispersion relation for the magnons can then be written as 

^=-4 7 (e( 1 +( ! t) 2 ))- (") 

This is the leading order behavior of the spectrum of the theory. There will in general be corrections to 
this formula involving higher powers of 4- However if we take the limit, where, -jV is small i.e the BMN 
limit 15 , then the above formula, which resulted from a one loop large N perturbative computation should 
be expected to capture the qualitative nature of the spectrum of the theory to leading order in J. 
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4.3 Bethe Ansatz for the Integrable SU(d') Sub-Sectors 

The large N one loop energy operator in the su(d') sector of the spin chains l|68ll can be written as l|79|l with 
a = — as the restriction to the su(d') sector simply amounts to dropping the trace operator from the 
spin chain Hamiltonian. The spins for this chain are in the fundamental representation of su(df). This spin 
chain is known in the literature as the generalized Heisenberg model. It is known to be integrable and the 
corresponding Bethe equations take on the form 

1^171,1 i 10L m .W \ TT ™'m,i ""m,j i ^^m-^m ^m,i ^"m',j < ^Otjn-Qm' (100) 

a m are the simple roots of the Lie algebra and w is the highest weight of the representation. The details of 
the derivation of these Bethe equations have been reproduced in the appendix. To use the Bethe equations 
and derive the results that put this generic case on the same footing as the previous results on the special 
case of the so(2) case let us analyze the case of the above equation corresponding to all the impurities of the 
1 type. In general, one can have as many types of impurities as the rank of the Lie algebra. However, it is 
enough to study the special case m = 1 to get a feel for the nature of the spectrum. In this case, callingu^ 
fli, the Bethe equations become 

' fij+i/2 \ J = yr flj - ^fc + i/2 

Hi-%/2) 11 Mj - - Mfc - i/2 
By writing the equation in terms of the momenta pi 

e~ tp = (102) 
(i + 1/2 

rather than the rapidities, we can write the Bethe equations as 

m 

J Pk J = TJ S (p k>Pj ), (103) 



e~' 



where 



1 -L. p«(Pl+P2) I OpiP2 

S( Pl ,p 2 ) = - .. ^ , - ■ (104) 



This is completely analogous to the equations for the so(2) spin chain, except that 6 = 1 in the present case. 
That is to be expected, because, the case of having impurities of the 1 type reduces the problem to the usual 
spin half Heisenbrg chain, which is the same as the XX Z chain with the value of the anisotrpy 5=1. As 
discussed in the appendix, the dispersion relation takes on the form 

m _^ 

E = Y^e(p i ) i e(p) = 2 C0 <P)- ( 105 ) 

As before, one can also take the limit of large J, in which case the magnons/impurities behave as free 
particles on a circle. In this limit the momenta are given by 

Pi = (106) 
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where rii are integers. They can take on any value, but must be constrained to obey the level matching 
or zero momentum constraint implied by the cyclicity of the trace i.e J2 i rn = 0. The dispersion relation, 
simplifies to 



If one allowed for other types of impurities, then the form of the Bethe equations would be more complicated, 
however, the qualitative facts about the nature of the spectrum given by the simplest case above do not 
change. The effect of the other types of impurities, in a chain of finite length would be to change the allowed 
values of the momenta. However, the dispersion relation would remain the same. Moreover, only u\ enters 
the dispersion relations, thus the only effect that the other types of impurities would have are indirect. 
Although the Bethe equations written out in their full form are rather complicated algebraic equations, the 
lessons learnt from the simplest case of 1 type of impurities carry over to the more general scenario as well. 

4.4 Enhancement of Integrability in the large J Limit 

From the previous discussion of the construction of integrable so(n) invariant spin chains with nearest 
neighbor interactions from the the standard 1Z matrix it follows that the form of the Hamiltonian is highly 
constrained from the requirements of integrability. Only Hamiltonians with very special relative coefficients 
between the permutation and the trace operator fulfill this requirement. Although, we were able to find 
rather large integrable subsectors for the so(d) spin chains, it is worthwhile to ask the question whether 
or not one can find integrable behavior in some limit of the spin chains without having to resort to the 
truncation to sub-sectors. If it is indeed possible, then we should be able to gauge the behavior of the 
complete spectrum of the theory albeit in some appoximation scheme. The simplification of the nature of 
the spectrum in the large J limits of the various exactly solvable sectors described in the previous sections 
suggests that the large J limit might be such an approximation. This is indeed the case. In this final section 
on the Bosonic matrix models of the paper, we shall formally demonstrate that the large J limits of the 
various so(d) spin chains studied so far can be mapped to integrable non-linear sigma models irrespective of 
whether the original spin chain is integrable or not. 

Since we are interested in spin chains with nearest neighbor interactions, and nearest neighbor interactions 
are nothing but the lattice Laplacian, the non-linear sigma models we expect to approximate the low lying 
spectrum of the spin chains will be free sigma models with an appropriate homogeneous space as its target 
space. However, most classical sigma models of this kind are known to be integrable. Thus, it is reasonable 
to expect that the low lying excitations of an so(d) invariant quantum spin chain with nearest neighbor 
interactions will be described by an integrable classical sigma model even if the original spin chain is not 
integrable. We shall now proceed to make this idea more precise. We begin by writing the 1Z matrix in a 
slightly different but equivalent form 



Pij,Kij are the operators that permute and trace over the spins belonging to the i and jth vector spaces 
respectively, h is the natural 'quantum' deformation parameter in the problem. An alternative way of 
writing the R matrix is as a d x d operator valued matrix. For example, TZij can be expressed as a matrix, 
each of whole elements is an operator on the i th vector space as, 




(107) 




(108) 



K(u)l 




— u 



1 



S tlv {i),n,u = l---d. 



(109) 
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In writing the R matrix in the second form, we have introduced the Weyl operators S^ <v (i), which act as 
\fi >< v\ on the i th vecotr space. The Weyl operators satisfy the standard commutation relations, 

[Sn V {i), S a p(j)] = K5ij {8 va S^ - 5^S av ) ■ (110) 

The relations satisfied by the transfer matrix, the so called RTT relations, are written most succinctly as, 

J2i, 2 (« - «)Ti(«)T 2 (t;) = T 2 {v)T x (u)Ri, 2 {u - v). (Ill) 

1,2 refer to two auxilliary spaces. It is more insightful to write this equation out in terms of the matrix 
elements of T. 

^[T ab (u),T cd (v)] + ^-^ (T ad {u)T cb (v) - T ad {v)T cb {u)) + 

— ^ r (8 a . b T lc (u)T ld {v) - S c , d T al (v)T bl (u)) = (112) 

rig — (u — v ) 

We shall now write the corresponding RTT relation in the large J limit. To do that we shall define 

n = j,u',i/ = j,j,s = j (ii3) 

In terms of the rescaled variables, the TZ matrix becomes 

n{u% v = 5„„J(i) + + ^J—S^ii) « 8^I(i) + + ^jS^ii) + O(-l). (114) 

J 2 

Needless to say that the S operators satisfy 

h' 

[S^(i),S af 3(j)] = — {SvaSpp - 5rfS av ) . (115) 

From the two preceding equations it is quite clear that the term in the R matrix that involved g' , and which 
was ultimately responsible for the g' dependence of the Hamiltonian, is lower order in j and that ^ plays 
the role of the deformation parameter in the large J limit. The RTT relations can be written out as 

hT ab (u'),T cd (v')] + -^— t (T ad (u')T cb (v') - T ad (v')T cb (u')) + 
u — v 

, A n ( S a.bTi c (u')T ld (v') - S c<d T al (v')T bl (u')) = (116) 

(U V ) 

We also observe that 



[T ab (u'),T cd (v')] = j{t ab (u'), t cd (v')} + O(-l) (117) 

where t, as a matrix is the same as T, except that its entries are to be thought of as ordinary functions and 
not operators. Thus to leading order in -j, the RTT relations may be replaced by 

{tabid), t cd {v')} + (t ad {u')t cb {v') - t ad {v')t cb {u')) + 

u — V 

7- (S a Mic(u% d {v') - 8 c4 t al {v')t b i{u')) = 0, (118) 



(«' - V 1 ) 
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This simply implies that the matrix elements of t satisfy the following Poisson brackets: 



{*(«') ? t(v')} = [r(u' - v'),t{u') ® (119) 

where, 

rfz-y) = (P-iO- (120) 

x-y 

To clarify the notation it is worth mentioning that equation l|119fl is to be thought of as expressing the 
equality of the actions of both sides of the equation on two copies of a d dimensional vector space [41]. The 
main implication of (|119(l is that Tr£ is the generator of an infinite family of conserved charges i.e 

{Trt(u),Trt(u)} = 0. (121) 

Our analysis of large J integrability of the continuum sigma models obtained from Ij68l69|l would be complete 
if we can show that the relevant sigma models are contained in the family of commuting charges that follow 
from the classical transfer matrix t. This is indeed the case. From a technical point of view, taking the 
continuum limit of the spin chain amounts to taking its coherent state expectation value in a long coherent 
state of length J and deriving an effective Hamiltonian in terms of the group parameters characterizing the 
coherent state. For so(d) valued spin chains with nearest neighbor interactions this has been worked out in 
detail injS], and we shall refer to that paper for the detailed descriptions of many of the results that will be 
used here. The coherent state construction of[S| when applied to (|68I69|I generate the following sigma model 
Hamiltonins 

H = a J dxTr (d x Md x M) . (122) 

a is the effective coupling of the sigma model. It is — s Iji for (|68[) and §^tjit for (|69[) . x is nothing but 
the contunuum coordinate along the length of the spin chain. M(x,t) is a d x d antisymmetric matrix and 
it satisfies the Poisson brackets 

{M(x) tJ1 M(x') kl } = 5(x - x')(8 jk M a (x) + 8 a M jk {x) - 6 tk M 3l (x) - 5 3l M lk (x)) (123) 

The matrix M is built out of group parameters that parameterize the so{d) coherent state and it can be 
written in terms of complex d vectors Zi that satisfy 

ZiZi = Z*Z* =0,Z*Z i = l (124) 

as 

Mij = ZiZ* -ZjZ*. (125) 

This parametrization along with the various various algebraic identities for M that it leads to have also been 
discussed in 9 . The equations of motions for M can be derived using the Poisson structure given above 

d t M = 4ad x [M, d x M] (126) 

The equations of motion are equivalent to the following flatness condition. 

[d t + At,d x + A x ] = 0, (127) 
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where 

A x = —M 

u 

A t = —[M,d x M]-^M (128) 
u u 

are the two components of the Lax connection. The Poisson brackets between the spatial components of the 
Lax connection can be written down as 

{A x (x,u), A x (x',u')} = S(x - x')[r(u - u'),A x {x,u) ® I + 1 ® A x {x' ,u')]. (129) 

The above equation expresses the Poisson bracket between the dynamical variables as a Lie bracket in 
so(d). This is the familiar Lie-Poisson structure characteristic of classical integrable two dimensional field 
theories . The immediate consequence of the Lie-Poisson structure is that 

t{u) = Pe$« AAx)dx (130) 

satisfies the Poisson brackets (I119|) . This of course implies the fact that Tr£(u) is a generating function of 
the integrals of motion for the Hamiltonian Ijl22(l . This establishes the integrability of the sigma model 
describing the large J limit of the spin chain. We thus see that all the nearest neighbor spin chains that 
arise as the one loop energy operators for matrix models with dynamically generated masses correspond to 
integrable sigma models in the large J limit. 

A cautionary remark is in order while discussing the continuum limits of the spin chains (|68l I69|) . The 
so(d) invariant spin chain Ij68|l is Ferromagnetic and the continuum limit described above does approximate 
the spectrum around the Ferromagnetic ground state. However, (|69|) is anti-ferromagnetic. Thus it is to 
be understood that the sigma model following from the above construction, when applied to the so(3) case 
approximates the spectrum of the spin chain around its highest energy state. The sigma model following 
from the continuum limit of the so(3) chain around the anti-Ferromagetic vacuum is a level two su(2) 
Wess-Zumino model, and we shall refer to |39lR*0"| for the details of that continuum limit. 



5 Mass Deformed Super-Membranes and Integrability: 

In the final section we shall analyze models of mass deformed supersymmetric quantum mechanics obtained 
from the dimensional reduction of minimally supersymmetric Yang-Mills theories in various dimensions. The 
resultant matrix models can be taken to be matrix regularized light cone descriptions of super-membranes 
propagating in plane wave type backgrounds in one higher dimension than the one in which the original 
gauge theory was formulated. The gauge theories in question are minimally supersymmetric Yang-Mills 
theories in spacetime dimensions six four and three. If one carried out a naive dimensional reduction of these 
Yang-Mills theories to one dimension, the resultant matrix quantum mechanics system would not have any 
quadratic mass terms. The matrix quantum mechanics models can nevertheless be interpreted as regularized 
versions of non-critical M(atrix) theory. This is very much along the lines of the original BFSS proposal for 
M(atrix) theory in eleven dimensions. The gauge theory in question there is of course AT = 1 SYM in D = 10. 
In the presence of enough supersymmetries, one cannot be confident of the success of the mechanism for 
dynamical mass generation that we applied to the case of Bosonic membranes earlier in the paper 6 . There 

6 In a recent paper 26 it has been argued that if one allowed for central extensions of the underlying SUSY algebra for the 
matrix models then mass gaps might be generated even in the presence of supersymmetry. It would of course be interesting 
to see if the methods of resummed perturbation theory can capture the mass gaps of systems posessing centrally extended 
supersymmetry 
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might be 'true' flat directions in the quantum mechanical system due to cancelations between Bosonic and 
Fermionic self-energy diagrams. 7 However, there does exist a second mechanism for adding mass terms in 
supersymmetric matrix models. One can add explicit mass terms to the matrix model Hamiltonian and 
ask for for a preservation of all the supersymmetries of the un-deformed model. This greatly constrains the 
permissible types of mass deformations, and all possible mass deformations of the matrix models obtained 
from dimensional reductions of minimally supersymmetric Yang-Mills theories in various dimensions were 
carried out in|21|. Sometimes, it is possible to interpret the massive matrix quantum mechanics as being 
related to the dimensional reduction of gauge theories, not on R m but on on R x M m_1 , where M is a m — 1 
dimensional compact manifold. For example the BMN matrix model can be regarded as the dimensional 
reduction of ftf = 4 SYM on R x S 3 . Obviously, the question as to whether or not there is always such 
an interpretation of massive supersymmetric matrix quantum mechanical systems possibly requires further 
study. Apart from the connection to Yang-Mills theories, such mass deformed matrix models can be regarded 
as light-cone Hamiltonians for super-membranes in plane wave like backgrounds, once again, in one dimension 
higher than the dimensionality of the original Yang-Mills theory |2T). Thus, they are quite naturally models 
for non-critical M(atrix) theory Hamiltonians. Thought of in another way, this is nothing but the extension 
of the relation of the BMN matrix model to a critical supermembrane theory^] in a plane wave background 
to non-critical dimensions. Taking a clue from the recently discovered integrable behavior in the BMN 
model [19] it is natural to ask, whether or not the matrix models for non-critical membrane correspond to 
integrable spin chains in the large N limit. That is the problem that we shall concern ourselves with in this 
part of the paper. As in the BMN model[T2|, a good starting point for the anlysis of integrability is provided 
by the Bosonic sectors of these models at one loop. Without integrabilty in the Bosonic sector at the one 
loop level, there is of course no hope of expecting integrability at higher loops. In what follows, we shall 
systematically analyze the Bosonic sectors of all models proposed by Kim and Park|22] at the one loop level 
at large N. 

The various mass deformed models obtained by Kim and Park|2lj can be classified into two broad 
catagories. 

I: Models that admit non-trivial supersymmetric configurations saturating the unitary bounds of the corre- 
sponding super algebras. 

II: Models that do not posses any non-trivial supersymmetric configurations 

In what is to follow, we first catalogue all the bosonic subsectors of the models presented inj2J|- We then go 
on to show that all the models that admit supersymmetric vacua also correspond to integrable spin chains 
at the one loop level. On the other hand, models that do not have non-trivial supersymmetric vacua, do not 
correspond to spin chains that are integrable, but they do possess integrable sub-sectors. 

5.1 J\f =8 Super Matrix Quantum Mechanics 

As mentioned before, the massive matrix model obtained from J\f = 1 SYM in D = 10 is nothing but the 
BMN matrix model and it has been studied in substantial detail in the recent past^Hj- The theory involves 
sixteen supercharges, and it shares many common features with Af — 4SYM on R x S 3 . Among the list of 
common features is the integrability of a rather large sub-sector of the matrix model, the so called su(2|3) 
sub-sector to the third order in large N perturbation theory[Hl]. At the one loop order, the complete theory, 

7 Even for supersymmetric matrix models, one can apply re-summed perturbation theory to estimate mass gaps at finite 
temperatures. Indeed, such an analysis for the eleven dimensional M(atrix) theory Hamiltonian has been carried out by Kabat 
and collaborators in 1251 . It might thus be possible to set up a systematic perturbation theory using a combination of spin-chin 
techniques and the mechanism of dynamical mass generation in the case of supersymmetric matrix models at finite temperatures 
as well. 
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without any truncation to specific sub-sectors also exhibits integrability 8 . The next lower dimension that 
admits a consistent formulation of Af = 1 SYM is six. The mass deformed matrix quantum models obtained 
from this Yang- Mills theory are of two types. The matrix model that has non-trivial susy vacuum solutions 
has as its Bosonic part 

H = Tr (hlilli + \{^) 2 X a X a + ^fX A X A - ±[X h Xtf + ip[X 3 , X 4 ]X^j (131) 

a = 3, 4, 5 and A = 1, 2, while i, j in the quartic interaction term run from 1 ■ • ■ 5 in the above equation. The 
mass deformation breaks the SO(5) symmetry of the original gauge theory down to 5*0(3) x SO(2). This 
models admits a supersymmetric configuration corresponding to circular motion in the (1,2) plane and a 
fuzzy sphere in the (4,5,6) directions. 

Xi = Rcos(-tn)X, X 2 = ifein(-i/i)T, X a = \^J a (132) 
6 6 3 

J are the standard so(3) generators satisfying 

[Ja,Jb] = itabcJc- (133) 

Clearly, states that involve only Bosonic excitations are eigenstates of the free part of the Hamiltonian, which 
is a sum of decoupled harmonic oscillators. It thus makes sense to compute the corrections to the tree level 
energies of the Bosonic excitations of this model. To do such a computation at the one loop level, one will 
have to take into account the contributions from the Fermionic part of the Hamiltonian, since the Fermions 
do run in loops. However, as far as the one-loop energy operators is concerned, the analysis of the BMN 
matrix models shows that contribution of the Fermionic terms to the one loop effective Hamiltonian, when 
written down in the spin chain language is proportional to the identity operator |191 122) . But since we are 
interested in measuring the energies with respect to a reference state, that contribution of the Fermions to the 
energy operator as we have defined it in the paper vanishes. This general observation is true of the present 
and all the other models that we study in this section. Since the symmetry algebra of the matrix model is 
so(3) x so(2), we can read off the contributions of these sectors to A 1 , using the basic relations 1)3 1132)1 . 

A -(2) = ^E(^M+i-^ + i)- (134) 

A ~<3) = ^2 E ( 2P W1 - 4i ^+l) ■ ( 135 ) 

The dimensional reduction of six dimensional super Yang-Mills theory also admits a second type of mass 
deformation, where the bosonic so(5) symmetry is broken to so(4) x u(l). The Bosonic Hamiltonian for this 
mass deformation is 



h = Tr^n l n l + i(|) 2 M + i(|) 2 ^Y A -i l x «, x tf) 

\i=l A=l a<b=l ) 



(136) 



8 An exhaustive account of the one-loop integrability of the BMN matrix model along with its relation to Af = 4SYM can 
be found inl231 
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In this mass deformation, there are no BPS solutions analogous to the fuzzy sphere solution of the so(3) x so(2) 
case. When one looks at A 1 for this model, the only non-trivial contribution to it comes from the so(4) part. 
As before, using (|31|) we can write down 

A S 1 o(4) = ^E(^M+i-^ + i). (137) 

The spins are now in the defining representation of so(4). 



5.2 jV =4 Super Matrix Quantum Mechanics 

Next in our catalog, are the matrix models obtained from a mass deformation of the dimensional reduction of 
minimal super Yang-Mills in D = 4. As in the previous case, one can do two consistent mass deformations |2l] . 
The Bosonic part of the Hamiltonian that admits a maximally supersymmetric fuzzy sphere configuration is 



H = Tr ( -11*11; + i(|) 2 Y a A a - I [A,, A,] 2 + ifi[X 1} X 2 ]X 3 ) ( l:>>N 



This matrix model corresponds to the so(3) symmetric part of the previous model and the corresponding 
fuzzy sphere vacua are also the same. 

X a = ^Ja, [Ja,Jb] = l^abcJc (139) 

The one loop energy operator for the scalars of this model is also precisely the same as the one obtained 
above. 

A -0) = ^2 E ( 2P M+i - 4A M+1 ) . (140) 

A second kind of mass deformation of the dimensional reduction of D = 4,J\f = 1 SYM results in a Bosonic 
Hamiltonian that has the so(3) symmetry broken down to so(2) x u(l). The Hamiltonian for the Bosonic 
part is 

H = Tr Qn,^ + ^ t i 2 (X 1 X 1 + X 2 X 2 + 4X 3 X 3 ) - i[X t , Xj} 2 ) (141) 

This model too admits a supersymmetric configuration preserving all of the four supercharges 21 . The 
supersymmetric directions correspond to 

Ai = Rcos(-tu)l, X 2 = Rsm(-tu)X, A 3 = (142) 
6 6 



A 1 -H 



E (Vu+i - p w) ■ ( 143 ) 



5.3 J\f = 2 Super Matrix Quantum Mechanics 

We now move on to the final case, which is that of dimensional reduction of three dimensional minimal super 
Yang-Mills. Unlike the previous cases, there is a unique mass deformed matrix model that one can obtain 
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from minimal super Yang-Mills in three spacetime dimensions. The Bosonic part of the Hamiltonian of this 
model is 



// = Tr ( in 4 n 4 + ^ 2 {x 1 x l + x 2 x 2 ) - \[x u x 3 ] 2 



(144) 



There is a maximally supersymmetric BPS configuration ins this case as well. It is given by the static 
solution 

Xi = Rcos(-tu)l, X 2 = Rsm(-tfi)I. (145) 
6 6 

Finally, we write down the one loop energy operator for the model, which takes on the form: 

A -(2)-^E(^M + i-^ + i). (146) 

Integrability: 

Having listed all the bosonic one-loop energy operators for the mass deformed models, we can just easily 
understand whether or not the spin chains are integrable or not. The basic results are summarized in the 
following table. 



Supercharges 



4 
4 
2 



Symmetry 
so(3) x so(2) 
so(4) x u(l) 

so(3) 
so(2) x u(l) 

so(2) 



Integrable Sectors 
Both so(3) and so(2) sectors 
su(2)(e so(4)) sector 
so(3) sector 
so(2) sector 
so(2) sector 



From the summary given above, we see that all the cases that do support non-trivial fuzzy sphere type of 
supersymmetric configurations give rise to integrable spin chains. The corresponding spin chains have so(3), 
so(2) or their product (corresponding to Af = 2, 4 and 8 respectively) as their symmetry groups. Integrability 
of the so(2) sectors of the Af= 8, 4 and2 cases follows directly from the realization of so(2) invariant nearest 
neighbor spin chains as xxz models discussed previously. We can simply use these results to derive the Bethe 
ansatz for the three so(2) cases of interest in the supersymmetric context. This simply amounts to replacing 
fi by # in the formulae (|96I97|) for the so(2) Bethe ansatz given previously. 

As far as the two so(3) invariant sectors of the Af = 8 and Af = 4 cases are concerned, they are integrable 
as well! From form of the most general integrable so(d) spin chain with nearest neighbor interactions given 
in l|78|) . we see that the so(3) chains given above have the correct form commensurate with integrability. 
To get the spectrum and the Bethe equations, one can use the known Bethe equations for the defining 
representation of so(3) EBI E3 EHj • Since so(3) has only one Cartan generator, the Bethe equations only 
allow for one type of impurity and they are very similar to the so(2) equations. When written out in terms 
of the rapidities they read as 

"■■'). (147) 



u k - i 



while the dispersion relation is 



e(u) 



36 1 
H 2 u 2 + 1 



(148) 



for the so(3) chains related to both the Af = 8 and Af = 4 cases discussed above. This provides the complete 
one-loop solution for the two so(3) sectors of interest. 
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The one example of a mass deformed model that does not have any non-trivial supersymmetric vacua 
corresponds to the case of the Af = 8 matrix quantum mechanics with symmetry group so(4) x The 
so(4) spin chain (|137|l is clearly not integrable. However, as in the case of the matrix models for Bosonic 
memebrane theories in flat spacetimes discussed previously in the paper, one can find integrable su(2) 
subsectors of the theory by looking at states formed by excitations corresponding to Z = X\ + iX^ and 
W = X2 + 1X4. The resultant spin chain is nothing but the spin one-half Heisenberg model, whose Bethe 
equations are 



This completes the discussion of the one- loop spectra of the Bosonic sectors of the models proposed in|21|. 
Whether or not, this understanding of integrability extends to the more general sectors of theme models for 
supermembranes and whether or not integrability is preserved at higher loops are of course open questions 
that might be of interest. 

Before concluding the segment on the analysis of the matrix models given in|21j. it is worth remarking 
that we have left out one particular example that was also presented in|21|. This last example corresponds 
to J\f = 1 matrix quantum mechanics and is related to the dimensional reduction of minimal SYM in D = 2. 
This matrix model, which has a single Bosonic and a single Fermionic degree of freedom cannot quite be 
studied using the formalism presented here, as one needs at least two Bosonic matrices to be able to make 
it work. However, a different approach, might be in order here. Supersymmetric matrix quantum mechanics 
systems, involving a single Bosonic matrix degree of freedom can often be mapped to supersymmetric analogs 
of the Calogero-Sutherland model, see for example 43 . It might be possible to utilize such a connection in 
this case leading possibly to some exact statements about this model. This possibility probably merits 
further investigation. 

6 Concluding Remarks 

In this final section, we briefly return to the analysis of Bosonic membranes and sketch out some directions 
for future investigations. In sections three and four we presented a mechanism for estimating the non- 
perturbative mass-gaps in the spectra of Bosonic membranes and also developed an expansion of the matrix 
model Hamiltonians around the quantum corrected effective potentials using the techniques of quantum spin 
chains. It is of course implied that we were expanding the Hamitonian around the trivial/oscillator vacuum 
around the effective potential. The low-lying excitations around the effective potential were shown to be 
well described by closed spin chains. In the D-brane picture, the closed spin chains could be thought of as 
the low-lying closed string like excitations of the membrane. Indeed, the sigma model for the spin chains 
derived as the continuum limits of the spin chains can be thought of as an effective string sigma model in a 
background provided by the membrane. The membrane point of view then provides us with an open-string 
interpretation of the sigma model. One could take this interpretation seriously and ask if there are analogues 
of bona-fide open string degrees of freedom in the membrane picture as well. This it indeed possible, if one 
looked at particular O(N) excitations of the matrix models. The closed spin chains arose in our analysis 
because we focussed on gauge invariant operators built out of traces. However, taking a cue from the analysis 




(149) 



The total energy is given by 




(150) 
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of the dilatation operator of J\f = 4 SYM, we could just as well have looked at states built out of determinants 
and sub-determinants|44| . Such operators can also be studied within the paradigm of quantum spin chains, 
except, one would have open spin chains i.e the analogs of the open string degrees of freedom to consider. 
Issues related to integrability of these spin chains tend to be rather subtle (see for example03]) and it might 
be of interest to investigate these degrees of freedom for bosonic membranes as well. 

The complete spectrum of a theory of membranes is of course much richer than what has been discussed 
in this paper. The richness of the spectrum has to do with the fact that the trivial/oscillator vacuum of 
the bosonic matrix models is only one of many other vacua that are also present in the theories that we 
considered. An example of a non-trivial vacuum in a matrix quantum mechanical system was already alluded 
to in the section on supersymmetric matrix models. These arc the so called fuzzy sphere vacua. The Bosonic 
membranes, coupled to Chern-Simons fluxes that were considered earlier in the paper also possess such 
vacua 9 . These vacua continue to be present around the quantum/mass corrected effective potentials as well. 
It is well known that the expansions of the quantum mechanical matrix models around fuzzy sphere vacua 
leads to gauge theories in 2 + 1 dimensions. For instance, D = 2 + 1, J\f = 4 SYM and its connection to the 
expansion of the BMN matrix model around a fuzzy sphere vacuum has been discussed at length in Hi ! 17| 10 . 
Put differently, quantum mechanical matrix models can be thought of as providing a regularized description 
of gauge theories in three dimensions. The regularization is accomplished by introducing non-commutativity 
in the spacial directions. Gauge theories in three dimensions of course have mass-gaps in their spectra. 
For example, the mass gap in the spectrum of pure Yang-Mils in three dimensions has been computed 
in |49l I5U1 1511 1521 I53| . It is thus plausible that the methods for estimating spectral gaps for matrix quantum 
mechanics that were presented in this paper can well be utilized to compute mass gaps and glueball masses 
for D = 2 + 1 Yang-Mills theories. Such a possibility was also pointed out in[3B] . The mass gap in D = 2 + 1 
Yang-Mills can be estimated in two different ways. In 49 , the gap was computed by applying re-summed 
perturbation theory to the gauge theory, while in |5UJ 1511 1521 IS*3| a strong coupling expansion was utilized 
within the Hamiltonian framework. The answers obtained from both the approaches are strikingly close to 
each other. To draw a parallel at a matrix model level, the analysis of[2Sj is very much along the lines of 
|50l 1511 l5*2l I53| . while, the re-summed perturbation theory used in this paper is closer in spirit to|49j. As 
mentioned previously, it would be indeed remarkable if these results can be applied to D = 2 + 1 Yang-Mills 
theories with or without supersymmetry. We hope to return to this problem in a future publication. 

Finally, it is worth mentioning, that membrane dynamics can also be used to compute certain glueball 
masses for gauge theories in spaces of finite volume|2I]- Since, we can go beyond the computation of the 
masses and organize the computation of higher excited states of Bosonic membrane theories as quantum 
spin chain computations, our results might be applicable to improve upon the glueball spectroscopy of gauge 
theories in spaces of finite volume. 
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9 For a detailed study of various vuzzy sphere vacua in matrix models of the type considered in the paper see |54| 
10 The connection between fuzzy /non-commutative spaces and gauge theories has a rather large literature devoted to it. For 
a pedagogical introduction we refer to !48| . 
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7 Appendix A: 



In this appendix wc outline the Bcthe equations for an su(n) invariant spin chain with Hamiltonian given 

by 

H=-J2Pu+i- (151) 
I 

This discussion is meant to provide the background material for various formulae related to Bcthe ansatzae 
throughout the paper. For a comprehensive discussion of su(n) Bethe ansatz techniques we shall refer to 
the original papers[3Sl |23 E21 EH] ■ 

The spins or the relevant chain are taken to be in the (n) representation of the algebra, and hence have 
n components. To find the spectrum of the Hamiltonian, one starts with the Lax operator that acts on 
7~ii ®V a . Hi is the 'one particle Hilbert space' associated with the site I and it is nothing but C n . V a is the 
auxiliary vector space, and as a vector space, it is also C n . The Lax operator is|38| 

Li, a = aWi, a + b{v)Pi, a (152) 

where 

a W = rrr, 6(/x) = ttt. (153) 

fi — 1/2 [j,— i/2 

One builds the transfer matrix from the Lax operator as 

Tj, a (n) = L Jta ---L ha . (154) 

This transfer matrix satisfies the Yang-Baxter algebra 

Tlabip. - v)Tj ja (p)Tjj,{v) = Tjfi{v)Tj^ a {[i)TZ ab (n - v). (155) 

The trace of the transfer matrix over the auxiliary space i(/i) = Tr a Tj >a (/i) is the generator of conserved 
commuting charges, i.e. 

[t{n),t{v)]=Q. (156) 
The conserved charges can be obtained in the standard fashion by expanding the transfer matrix around 

° l=l {£) M^WOriU-,. (157) 

In particular, 

H = O 1 — J. (158) 
The momentum is given, in terms of the Bethe roots, by 

M + i/ 2 ' 



= ^ (159) 



while the dispersion relation 

e(p) = -2cos(p) (160) 

becomes 
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The Bethe equations are 



'hn. i "~ ^ ^m,i ^m,j \ if™Tw™m ""m,i Um' ,j > ^^m-^m' (162) 

. y . ^mA *^m*^7n , / / - ^m.i ^m\j ^Q! m .Q; m / 

j^ti J m ^m j^i J 

a m are the simple roots of the Lie algebra and w is the highest weight of the representation. The eigenvalues 
of the Hamiltonian in terms of the solution of the Bethe equations requires further details of the Bethe ansatz, 
and we shall simply quote the result here while referring to 35, 36]|3ZIEH] for a more complete derivation. 
The total energy, is given in terms of the dispersion relation by 

III 

E = -J2<u ld ). (163) 
j=i 

It is important to note however that only roots of the first type enter the dispersion relation, while the other 
roots only affect the energies indirectly. 

To be able to use the Bethe equations, one also needs information about the roots and weights for the 
Lie algebra, which we also list below for the sake of completeness. For su(n), the simple roots are given by 

a m = v m - v m+ i, m = 1 • • • n — 1, (164) 

while the n — 1 dimensional weights are given by 



The weights satisfy 



while 



[uj] m = = ( 2^ 5 jtk - mS j>m+1 ] . (165) 

■w 2TTl\J 



-^+o<W', (166) 
In I 



5mm 1 — ^^m.m' + l- (167) 



The heighest weight in the defining representation w = V\. 
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